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Abstract. 

Earth-like planets have anelastic mantles, whereas giant planets may have anelastic cores. As for the fluid parts of 
a body, the tidal dissipation of such solid regions, gravitationally perturbed by a companion body, highly depends on its 
internal friction, and thus on its internal structure. Therefore, modelling this kind of interaction presents a high interest to 
provide constraints on planet interiors, whose properties are still quite uncertain. 

Here, we examine the equilibrium tide in the solid central region of a planet, taking into account the presence of a 
fluid envelope. We first present the equations governing the problem, and show how to obtain the different Love numbers 
that describe its deformation. We discuss how the quality factor Q depends on the rheological parameters, and the size of 
the core. 

Taking plausible values for the anelastic parameters, and examinig the frequency-dependence of the solid dissipation, 
' we show how this mechanism may compete with the dissipation in fluid layers, when applied to Jupiter- and Saturn-like 

planets. We also discuss the case of the icy giants Uranus and Neptune. 
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1 Introduction 

Once a planetary system is formed, its dynamical evolution is governed by gravitational interactions between its com- 
ponents, be it a star-planet or planet- satellite interaction. By converting kinetic energy into heat, the tides pertub their 
orbital and rotational properties, and the rate at which the system evolves depends on the physical properties of tidal 
dissipation. Therefore, to understand the past history and predict the fate of a binary system, one has to identify the 
dissipative processes that achieve this conversion of energy. Planetary systems display a large diversity of pl anets, with 
telluric planets having anelastic mantles and giant planets with possible anelastic cores (lUdry & Santosll2007l) . Since the 
tidal dissipation is closely related with the internal structure, one has to investigate its effects on each kind of materials 
that may com pose a planet. Studies have been carried out on tidal effects i n fluid bodies such as stars and envelopes of 
giant planets (lOgilvie & Linll2004l l20Q7b fo gilvid [20Q9L iRemus et al.lEoi2h . However, the planetary solid regions, such 
as the ma ntles of Earth-like planets or the ro cky cores of giant planets may also contribute to tidal dissipation (see for 
example fefroimskyll2012tlR^mus et al]|2012 ). We explore here the tidal dissipation in these solid parts of planets. 



2 The system 

Two-layer model. - We will consider as a model a two-bodies system where the component A, rotating at the angular 
velocity Q, has a viscoelastic core of shear modulus yu, made of ice or rock, surrounded by a fluid envelope, such as an 
ocean, streching out from core's surface (of mean radius R c ) up to planet's surface (of mean radius R p ). Both core and 
envelope are considered homogeneous, with constant density p c and p respectively. This model is represented on the left 
panel of Fig.[T] 
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Configuration. - We undertake to describe the tide exerted by B (of mass m#) on the solid core of A, when moving 
in an elliptic orbit around A, with eccentricity e, at the mean motion oj. Since no assumption is made on the B's orbit, 
we need to define an inclination angle / to determine the position of the orbital spin of B with respect to the total angular 
momentum of the system (in the direction of Z R ) wich defines an inertial reference plane (Xr, Y r ), perpendicular to it. 
The spin axis of A then presents an obliquity s with respect to Z R . Refer to the right panel of Fig. [T]for a synthetic 
representation of the system configuration. 



A(M A ) 





Fig. 1. Left: the system is composed by a two-layer main component A, with an homogeneous and incompressible solid core and 
an homogeneous static fluid envelope, and a point-mass perturber B orbiting around A. Right: B is supposed to move on an elliptical 
orbit, inclined with respect to the inertial reference plane (X R , Y R ). The equatorial plane of A (X E , Y E ) is also inclined with respect to 
this same reference plane. 



To treat the complexity of the two-layer problem, we follow the methodology of Dermott dl979b , 
3 Tidal dissipation of the core in the case of a two-layer planet 

Definition. - The tidal perturbation exerted by B on the solid core of A results on one hand in its deformation, and 
on the other hand in the dissipation of the tidal energy into heat leading to a lag angle 6 between the line of centers and 
the tidal bulge. This process can be modeled by the complex secon d-order Lo ve numb er k 2 defined as t he ratio of the 
erturbed gravific response potential ®' over the tidal potential U dBiotlll954l see also lTobiell2003l and iHenning et al.1 



pe 
2C 



2009). Its real part represents then the purely elastic deformation of the potential of the core ($') while its imaginary part 
accounts for its anelastic tidal dissipation. 

In practical calculations, we first have to develop U (and therefore ®') on spherical harmonics (F™), each term having 
a wide range of tidal frequencies (T2,m,p,q = (2 - 2p + q) qj - m£i, for ( m, p, q) e [f - 2, 2] x [f O, 2H x Z, resulting from the 
expansion of U on the Keplerian elements using the Kaula transform dKaulal ll962. see also Ma this & Le Poncin-Lafittd 
12009b , Thus, the complex Love number &2 depends on the tidal frequen cy and the rheology of the core, and so does the 
quality factor Q which quantifies the tidal dissipation (see for example lTobidl2003l) 

Q (v,o~ 2 , m , p , q ) = — 7=- ^f- , where k 2 (fi, o- 2 ,m, P , q ) = — ( = M e l l 2d ^, m , P , q )\ ? (3>1) 

CT2,m,p,q)\ U KO~2,m,p,q) 

where the quantity jl = fii +ifi2 = 9 l9fi P is the complex effective shear modulus, linked with the anelasticity (and thus 

^ Pc 8c 

the rheology) of the planet's core and its gravity g c . 

Case of a two-layer planet. - Acting as an overload on the solid core, the fluid shell, deformed by the tide, modifies 
both the tidal deformation and dissipation of the core. The second order Love number k 2 takes then a different form than 
in the fully- solid case 



0~2,m,p,q) = 



(5 + ^l C+ ^) + ^ 



- iA Dp2f , (3.2) 



where a, A, B, C and D account for the planet's internal structure through the ratios of radii §r and densities — . 

K p Pc 
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9 ^2(^2,m,p,q) 2 

4a 2 A 2 D 2 



1 + 



(B + jJL X {o-2,m,p,q)) + !*1 (^m,^)) 



V2(CT2,m,p,q) 2 



(3.3) 



depends on the core's parameters (its size, density and rheological parameters) and the tidal frequency. Moreover, to 
derive this expression of Q, no assumption has been made on the rheology of the core, except that it is linear under the 
small tidal perturbations (i.e. core's material obeys the Hooke's law). Hence, it is valid for any linear rheological model. 



Comparison with observations. - To confront our model with observations, we need to introduce the global dissipation 
factor, corresponding to a reseating of the previous one to the planet surface and thus involving the second-order Love 
number at the surface of the planet 



k 2 (R P ) 



k 2 (R c ) 



xQ. 



(3.4) 



Moreover, we need to choose a model to represent the way the core's material responds to the tidal perturbati on, i.e. 
a rheo logical model. Thus, from now on, we assume that the core behaves like a Maxwell body (see, for example jTobie 
20031) . 



4 Application to giant planets 



Application to gas giants. - Using astrometric data covering more than a century, lLainey et al.l (120091 120121) suc- 
ceeded in determin ing from observations the tidal dissipation in Jupiter and Saturn: namely, <2j up iter = (3.56 ± 0.56) x 10 4 
( Lainev et ail 2009). and Qsatum = (1.682 ± 0.540) x 10 3 determi ned bvlLainev et al.l d2012l) . Note that such high dissipa- 
tion is required by the formation scenario of Saturn's system of ICharnoz et al. ( 201 lb , in which the mid-sized satellites 
are formed at the edge of the rings. These values, which seem to be in agreement with other observations related to 
Jupiter's and Saturn's systems (see the corresponding references cited just above), are lower of up to one order of mag - 
nitude than what was expected by previous formation scenarios (see, for example, bfoder & Pealdll9 81: Sinclair 1 983 ). 
and even lowe r than what the most up-to-date models of fluid tidal dissipation predict (see, for example, lOgilvie & Linl 
2004; Wu 2005). Then, the question arises on the role of the poss ible solid central regions as sources of dissipation. Since 
the composition of giant planets cores is poorly constrained ( GuillotlE005h . we explore in Fig. [2] the tidal dissipation of 
Jupiter's and Saturn's core for a large range of values of the viscoelastic parameters considering the Maxwell rheological 
model. The other parameters (planet and core sizes and masses) are indicated in the legend. 
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Fig. 2. Dissipation quality factor Q eS as a function of the viscoelastic parameters G and 77, of a two-layer gas giant, using the Maxwell 
model. Left: for a Jupiter-like planet at the tidal frequency of Io. Right: for a Saturn-like planet at the tidal frequency of Enceladus. 
The red dashe d line indicates the value of geff = {(3.56 + 0.56) x 10 4 , (1.682 + 0.540) x 10 3 } (for Jupiter and Saturn, respectively) 
determined bv lLainev et all (2009. 2012). The blue lines corresponds to the lower and upper limits of the reference values taken by the 
viscoelastic parameters G and rj for an unknown mixture of ice and silicates. We assume the values of R p = {10.97, 9.14} (in units of 
R @ X M p = {317.8,95.16} (in units of M e ),R c = {0.15, 0.26} x R p , and M c = {6.41, 18.65} (in units of M©). 
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In l2004l studied tidal dissipation in rotating giant planets resulting from the excitation by the tidal 

potential of inertial waves in the convective region. Taking into account the presence of a solid core as a boundary 
condition for the reflexion of inertial waves, they obtained a quality factor geff ~ 5 x 10 5 . 

The present two-layer model proposes an alternative process that may reach the values observed in lLainev et all [2009. 



20121) . depending on the viscosity 7] and the stiffness G. 



To explain the tidal dissipation observed in the gas giant planets of our Solar System, all processes have to be taken 
into account. 



Application to ice giants. - As in gas giants, the standa rd three-layer mode l s for the interior structure of ice g iants 
predict the presence of a solid rocky core (see, for example, Hub bard et al1ll99ll:IPodolak et alJ lT995: Guillot 1999). But 
it still remains an incertitude on the phase state of the intermediate "icy" layer located between the rocky core and the 
convective atmosphere. Considering recent three-dimensi onal simulations of Neptune's an d Uranus' dynamo tha t predict 
that this region is a stably stratified conductive fluid one ( Stanley & Bloxhamll2004 120061) . iRedmer et alJ ( 201 lb studied 
the electric conductivity of warm dense water taking into account the phase diagram of water. Their results infer that part 
of this shell is in the superionic state, i.e. a two-component system of both a conducting proton fluid and a crystalline 
oxygen solid, and extends to about 0.42-0.56 of the planet radius. Thus, it seems reasonable to assume for our two-layer 
model that the solid central region extends from the rocky core surface up to somewhere in the superionic shell. 

We explore in Fig. [3]the tidal dissipation of Uranus' and Neptune's core for a large range of values of the viscoelastic 
parameters, considering the Maxwell rheological model, for different core sizes. 
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Fig. 3. Dissipation quality factor geff as a function of the viscoelastic parameters G and 77, of a two-layer ice giant, using the Maxwell 
model. Top: for a Uranus-like planet at the tidal frequency of Miranda, with three different core sizes R c = {0.12, 0.22, 0.32} x R p . Bot- 
tom: for a Neptune-like planet at the tidal frequency of Triton, with three different core sizes R c = {0.14, 0.26, 0.32} x R p . The red and 
orange das hed lines indicate, respectively, the lowest and highe st values of Q eff from formation scenarios: Q eS = {5 x 10 3 ,7.2 x 10 4 } 
for Uranus fcavrilov & Zharkovlll977l : lGoldreich & SoteJl96ot) and g eff = {9 x 10 3 , 3.3 x 10 5 } for Neptune dZhang & HamiltorJl2008l : 
Banfi eld & Mur ray 1992). The yellow dashed line indicates the value of Q eff from a study of Neptune's internal heat: Q eS = 1.7 x 10 
(Trafton 1974). The blue lines corresponds to the lower and upper limits of the reference values taken by the viscoelastic parameters G 
and 77 for an unknown mixture of ice and silicates. We assume the values of R p = {3.98, 3.87} (in units of /?©) and M p = {14.24, 16.73} 
(in units of M©). The core mass is obtained by integration of the density profiles o f Helled etdl bOllh up to a given core size. 
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5 Dynamical evolution 

Due to dissipation, the tidal torque has non-zero average over the orbit, and it induces an exchange of angular momentum 
between each component and the orbital motion. This exchange governs the evolution of the semi-major axis a, the 
eccentricity e of the orbit, th e inclination I of the orbital plane , of the obliquity s and that of the angular velocity of each 
component (see for example Mathis & Le Poncin-L afitte 2009). Depending on the initial conditions and on the planet/star 
mass ratio, the system evolves either to a stable state of minimum energy (where all spins are aligned, the orbits are 
circular and the rotation of each body is synchronized with the orbital motion) or the planet tends to spiral into the parent 
star. 



6 Conclusion 



Our evaluations reveal a much higher dissipation in the solid cores of planets than that found by|5 gilvie & Lin] ([2004) 
for the fluid envelope of a planet having a small solid core. These r esults seem to be in good agreement with observed 
properties of Jupiter's and Saturn's system dLainev et a l. 2009, 20121) . In the case of the ice giants Uranus and Neptune, 
too much uncertainties remain on internal structure to give an order of magnitude, other than a minimum value, of tidal 
dissipation in the solid regions, which constitutes a first step in the tudy of such planets. 



This work was supported in part by the Programme National de Planetologie (CNRS/INSU), the EMERGENCE-UPMC project EME0911, and the 
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